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Abstract 

Charge exchange spin-dipole (SD) excitations of 90 Zr and 208 Pb are studied by using a Skyrme 
c5 ' Hartree-Fock(HF) + Random Phase approximation (RPA). The calculated spin-dipole strength 

distributions are compared with experimental data obtained by 90 Zr (p,n) 90 Nb and 90 Zr (n,p) 
90 jsj]-, reactions. The model- independent SD sum rule values of various Skyrme interactions are 
studied in comparison with the experimental values in order to determine the neutron skin thickness 
of 90 Zr. The pressure of the neutron matter equation of state (EOS) and the nuclear matter 
symmetry energy are discussed in terms of the neutron skin thickness and peak energies of SD 
strength distributions. 

PACS numbers: 24.30. Cz, 21.10.Pc, 23.20.Lv, 21.60.Jz 
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I. INTRODUCTION 



The relationship between the neutron matter equation of state (EOS) and the neutron 
skin thickness has been studied extensively by using the Skyrme Hartree-Fock (HF) model, 
a relativistic mean field (RMF) model l[ [2], [J. The neutron matter EOS is essential for 
studying the properties of neutron stars, e.g., their size j^J. It is also known that isovec- 
tor nuclear matter properties, including the symmetry energy, correlate strongly with the 
neutron skin thickness in heavy nuclei 

m. 

Elastic electron scattering has provided accurate data on the charge distributions of 
nuclei. Several experimental attempts have been made to measure neutron distributions, 
for example, by proton elastic scattering [?], s|, Q and by inelastic alpha scattering to 
giant dipole resonance excitations . However, empirical results of neutron skin thickness 
obtained by proton scattering are controversial and do not agree with each other even within 
experimental error. The accuracy of empirical data on neutron distributions from giant 
resonance experiments is also rather poor, insufficient to extract accurate information on 
the neutron matter EOS. One promising tool for studying neutron distributions is the parity 
violation electron scattering experiment [ijj]. Unfortunately, no data on parity violation 
electron scattering experiments is available so far. 

The mo del- independent sum rule strength of charge exchange SD excitation is directly 
related to information on the neutron skin thickness 13[] . Recently, SD excitations were 



studied in 90 Zr by the charge exchange reactions 90 Zr(p,n) 90 Nb [14[ and 90 Zr(n,p) 90 Y 15], 

and the model-independent sum rule strengths for the SD excitations were extracted in Ref. 

[r| by using multipole decomposition (MD) analysis 17]. The charge exchange reactions 

( 3 He,t) on Sn isotopes were also studied to extract the neutron skin thickness 18| . However, 

one needs the counter experiment (t, 3 He) or (n,p) on Sn isotopes in order to extract the 

model-independent sum rule value from experimental data. This counter experiment is 

missing in the case of Sn isotopes. 

It is known that the SD strength has almost the same amount of contributions to neutrino 

PI 

reactions as that of the Gamow- Teller strength [19]. The Pb target is considered to be 
the most promising candidate for detecting the heavy-flavor neutrinos from the supernovae. 
Thus, it is quite important to study the SD strength in the Pb target for a precise evaluation 
of the cross-sections of charge-induced neutrino reactions. 
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In this paper, we study the SD excitations and the neutron skin thickness by using the HF 
and HF+ random phase approximation (RPA) with Skyrme interactions. As a theoretical 
model, the HF+RPA model has been extensively applied to giant resonances in a broad 
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region of mass table 
charge exchange excitations 
GT and SD states in 48 Sc and 90 Nb 
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he same model was used for the study of spin-dependent 
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26|. It was shown that the model successfully predicts 
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251 ] . First, we calculate the SD states in nuclei 



with mass A=90 and 208 by using the charge exchange HF + RPA model with various 
Skyrme interactions. We will compare calculated results of SD strength distributions with 
empirical data obtained by charge exchange (p,n) and (n,p) reactions on 90 Zr. The sum 
rule values are also compared with the empirical values in Section 2. Next, the correlations 
between the neutron matter EOS and the SD sum rules are studied in the Skyrme HF 
model. We will discuss the neutron matter EOS by using the experimental SD data and 
other empirical information on the neutron skin. This paper is organized as follows. In 
Section 2, the SD strength of the HF+RPA calculations is presented for both the £_ and t + 
isospin channels on 90 Zr and 208 Pb. The calculated results are compared with experimental 
results of 90 Zr(p,n) 90 Nb and 90 Zr(n,p) 90 Y reactions. We study the correlations between the 
sum rules of SD strength and the pressure of neutron matter EOS in Section 3. A summary 
is given in section 4. 



II. HF+RPA CALCULATIONS OF SD STRENGTH 

The operators for SD transitions are defined as 

S± = £4<£r<l?(ft) (1) 

with the isospin operators £3 = t z ,t± = (t x ± it y ). The model-independent sum rule for the 
A— pole SD operator S± = J2i^± r «[°" x ^i(^)] A can be obtained as 

S^-S* = ^ I (i I I 0) | 2 - I I S+ I °> I' 

i&all idall 

= (0 I [S x _, S x + ] I 0) = ^^W 2 >n - Z(r\). (2) 



The sum rule for the spin-dipole operator ([I]) then becomes 
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S--S + = - S + ) = Un{tX - Z{r\). (3) 



It should be noted that the sum rule is directly related to the difference between the 
mean square radius of neutrons and protons with the weight of neutron and proton numbers. 
We adopt four Skyrme interactions, namely, SIII, SGII, SkI3 and SLy4, for the HF+RPA 
calculations. The Landau parameters and nuclear matter properties of these interactions are 
shown in Table [H For the spin-isospin excitations, the value G' plays the important role of 
determining the collective properties of the excitation [22[|. The RPA equation is solved using 
the basis expanded by the harmonic oscillator wave functions up to the maximum major 
quantum number of N max =10 for 90 Zr and N max =12 for 208 Pb. The HF calculations are 
performed without the spin-gradient terms (J 2 terms) since the adopted Skyrme interactions 



have been fitted without them 



22j, |23|], but the RPA calculations incorporate the spin- 



gradient terms. The two-body spin-orbit and two-body Coulomb interactions are neglected 
in the RPA calculations. We also performed the continuum HF+RPA calculations with one 
of the interactions and found essentially the same strength distributions as in the present 
calculations except the width due to the coupling to the continuum 21] . The calculated 
results are smoothed out by using a weighting function, p: 

dB(SD) ave f dB(SD) 
dE x 

where the weighting function is defined as 



/ d -?£p-P(K - E r )dK (4) 



P(K - E.) = ^ (£ ,_ £ j 2 /2 +(A/2)2 (5) 

taking the width parameter A. In the present calculations with the discrete basis, the SD 
strength is given by 

d -^l = J2B(SD;E t )5( El -E' x ). 



A. Charge exchange SD excitations of J0 Zr 

The HF calculations are performed by using four Skyrme interactions in Table [H The 
proton, charge and neutron radii of 90 Zr are listed in Table [TT] together with the sum rule 
values AS = S_ — S + calculated through the analytic equation flSj). By using the same HF 
wave functions, the charge exchange RPA calculations give the SD strengths in 90 Nb and 90 Y 
excited by the t±raY\{f) operators from the parent nucleus 90 Zr, as shown in Figs. [1] and 
[2j The experimentally obtained distributions of the SD strengths are also plotted in Figs. 
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TABLE I: Landau parameters, effective mass m* and symmetry energy J of Skyrme interactions 





SIII 


SGII 


SkI3 


SLv4 


^0 


0.309 


-0.235 


-0.318 


-0.273 


U 


0.862 


0.733 


0.653 


0.818 


Gq 


0.052 


0.014 


0.569 


1.120 




0.457 


0.509 


0.203 


-0.138 


Fi 


-0.709 


-0.646 


-1.269 


-0.926 


1 


0.490 


0.521 


-0.843 


-0.399 


Gi 


0.490 


0.612 


1.33 


0.279 


G[ 


0.490 


0.432 


0.65 


1.047 


m* jm 


0.76 


0.78 


0.58 


0.69 


J(MeV) 


28.1 


26.9 


34.8 


32.3 



[U and [2j The experimental SD strength distributions for the t_ and the t + channels were 
obtained from the 90 Zr(p,n) 90 Nb and the 90 Zr(n,p) 90 Y data, respectively, by performing MD 
analysis 16]. A comprehensive description of the MD analysis can be found in ref. [1 71 ] . 



TABLE II: Proton, neutron and charge radii of 90 Zr. The charge radius is obtained by folding the 
proton finite size. The sum rule values AS" = S~ — S+ of spin-dipole excitations are calculated 
by Eq. ([3]) with the HF neutron and proton mean square radii. Experimental data on the charge 



radius are taken from ref. 



271 ] . The experimental values r n — r p are taken from 



m. 



The radii 



are given in units of fm, while the SD sum rules are given in units of fm 2 . 





SIII 


SGII 


SkI3 


SLy4 


exp 


r p 


4.257 


4.198 


4.174 


4.225 




T c 


4.321 


4.263 


4.240 


4.290 


4.258±0.008 


T n 


4.312 


4.253 


4.280 


4.287 






0.055 


0.055 


0.106 


0.064 


0.09±0.07 [7], 0.07±0.04 [16J 


AS 


146.7 


142.9 


156.9 


146.9 





In general, the t_ SD strength distributions for 0~ and 1~ states in 90 Nb are concentrated 
in one state at E x ~30MeV, having a large portion of the non— energy weighted sum rule 
(NEWSR) strength, while those for the 2~ states are separated into two dominant peaks, as 





FIG. 1: (Color online) Charge exchange SD strengths for the operators = ^i^- r i\ a x ^i(n)]^ 
calculated by the HF+RPA model with the Skyrme interactions (a) SIII, (b) SGII, (c) SkI3 and (d) 
SLy4. The excitation energy is referred to the ground state of the parent nucleus 90 Zr. The dotted, 
dashed and long— dashed lines show the SD strengths of A = 0~, 1~ and 2 _ , respectively, while 
the solid curve shows the sum of three multipoles. The SD strength is averaged by the weighting 
function ([5]) with the width A=lMeV. The experimental data shown by the black dots are taken 
from ref. [la ]. 

shown in Fig. [IJ The 0~ peak appears at E x ~30MeV, having 73%, 65%, and 58% of the 
NEWSR value for the SIII, SGII and SkI3 interactions, respectively. The calculated results 
for 1" states show a peak at E x ~29MeV having 50%, 59% and 48% of the NEWSR value 
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£,(MeV) 



FIG. 2: (Color online) Charge exchange SD strengths for the operators = r A a x ^i(n)]^ 

calculated by the HF+RPA model using the Skyrme interactions (a) SIII, (b) SGII, (c) SkI3 and 
(d)SLy4. The excitation energy is referred to the ground state of the parent nucleus 90 Zr. The SD 
strength is averaged by the weighting function (0) with a width of A=lMeV. The experimental 
data shown by the black dots are taken from ref. [la |. See the captions to Fig. [TJfor details. 

for the SIII, SGII and SkI3 interactions, respectively. The three results in Figs, Q] (a), (b) 
and (c) show the 0~ peak at a very similar excitation energy, while the values of NEWSR 
are somewhat different. The same is true for the 1 ~ peak in the three results. For the SLy4 
interaction in Fig. QJd), the 0~ and \~ peaks appear at about 3 MeV lower than the other 
three results, having 76% and 68% of the NEWSR, respectively. This is due to the negative 
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TABLE III: Peak energies and the average energies of charge exchange SD excitations in the A=90 
nuclei obtained by the self-consistent HF+RPA calculations: t_ in 90 Nb and t + in 90 Y. The average 
energy is calculated by the ratio of EWSR to NEWSR: E(MeV)=mi/m . See the text for details. 





t- 


t+ 




E peak {MeV) 


E(MeV) 


E peak (MeV) 


E(MeV) 


SIII 


28.5 


25.7 


13.5 


10.9 


SGII 


27.7 


26.7 


11.7 


9.47 


SkI3 


29.3 


28.2 


12.8 


11.6 


SLy4 


26.1 


24.9 


11.4 


10.5 



value of the Landau parameter G' in SLy4 for the spin-isospin channel. The dominant 
configurations of the collective 0" and 1~ states are the (Trl/ig^lgg^) and (7ilg 7 / 2 vlfy 1 2 ) 
configurations. For the 2~ excitations, the number of p-h configurations is larger than those 
of 0~ and and therefore the strength is fragmented in a wider energy range compared 
with 0~ and 1~ excitations. There is a small low-lying peak with J n = 2" at E x = 12.4 
(14.1) MeV with 10.0 (9.0)% of the NEWSR value in the case of the SIII (SGII) interaction. 
This state is mainly due to the 7iT #9/2^1/^2 configuration. The major strengths are found 
in the two peaks around 21 and 27 MeV in both the SIII and SGII results. The strength 
around E x = 21 MeV exhausts 50(41)% of the NEWSR value, while the peak around E x = 
27 MeV exhausts 30(37)% of the NEWSR value for the SIII (SGII) interaction. The peak 
energies in the two results are similar, while more SD strength is shifted to the peak around 
E x — 21 MeV in the case of the SIII interaction. The main configurations of the higher peak 
at E x = 27 MeV are the same as those of the 0~ and 1~ peaks, namely, (nlhg^vlgg^) and 
(771^7/2^1/^2)- On the other hand, the main configurations of the peak around E x = 21 
MeV are (7rl/in/ 2 ^l5' 9 ~/ 2 ), (7r2d5/ 2 ^2pr, 1 2 ) and (7r2rf5/ 2 z/2p~ // 1 2 ). The 2 _ strength distributions 
of the SkI3 and SLy4 interactions are somewhat different than those of the SIII and SGII 
interactions. There is no isolated low energy peak in the results for the SkI3 and SLy4 
interactions. Three large peaks are seen at E x = 23.5, 26.5 and 29.5 MeV together with 
several small peaks, while the two peaks at 20 and 26 MeV exhaust most of the strengths 
in the case of SLy4. 

The SD strengths calculated by the SD operator = r i\ a x ^i(^i)] A are shown in 
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Fig. [2](a), (b), (c) and (d) for the SIII, SGII, SkI3 and SLy4 interactions, respectively. The 
strength distributions are divided into two energy regions: a broad bump below 10 MeV 
and a peak around E x = 13 MeV. The strengths below 10 MeV are due to the \ n =l~ and 
2~ states, while the high energy peak is induced mainly by the 1~ states. The large 0~ 
strength is also found just above the high energy 1~ peak. The summed NEWSR values 
of all multipoles below 10 MeV are almost equal to the strength of the high energy peak 
around E x = 13 MeV in the case of the SIII and SGII interactions. The high energy peak 
of the SIII interaction in Fig. [2(a) is about 2 MeV higher than those of SGII and SLy4, as 
listed in Table HTT1 The main configuration for the high energy peaks with A 7r =0 _ l - and 
2~ is the (ulg 7 /2^fy 2 ) excitation. For the low energy 1~ strength, the (b , 2d 3 / 2 'n'2p~^ 2 ) and 
(^1^7/2^1/^2) configurations play the dominant roles. The (b / 2d 5 / 2 7i2p~^ 2 ), (b / 2d 5 /27i2p~^ 2 ) 
and (z/3si/2vr2p3 // 1 2 ) configurations have a large contribution in the low energy peak of A 71 " = 
2~. The (^2^3/2^1/^2) configuration contributes substantially to the high energy 2~ peak 
together with the (ulg 7 /27ilfy 1 2 ) configuration. The large spread in the distributions of the 
SD strengths in Figs. [Hand [2] are due to the fact that the p — h excitations are very different 
in unperturbed energy. Thus, the collision-less Landau damping effect plays an important 
role in the large observed width of SD resonance, while the coupling to the continuum 
plays a minor role. The coupling to 2-particle-2-hole (2p-2h) states were shown to increase 



28|]. 



substantially the width of the main peak of the £_ SD excitations of 90 Zr in ref. 

The energies of the main peaks E pea k are tabulated in Table UTTl along with the average ex- 
citation energies, which are calculated by the ratio of the energy-weighted sum rule (EWSR) 
mi to the non-energy weighted sum rule (NEWSR) mo, E= mi /mo. The E is always lower 
than the E pea k because of the low energy peak in the excitation spectra. For the t_ response, 
the SkI3 interaction gives the highest excitation energy for the peak, while the SLy4 is the 
lowest. Notice that the energy of SkI3 is the highest due to the small effective mass m*/m, 
while the negative Landau parameter G' is responsible for the fact that SLy4 yields the 
lowest energy value in Table IIII1 The general trend of the average excitation energy E is the 
same for the t + response. The SIII, however, gives a somewhat higher energy for the E pea k 
than SkI3 does. 

The calculated results of SD strength are shown in Fig.[3]together with the experimentally 



obtained distributions of the SD strengths 



161 ] . The spectra for the t + channel are shifted 



by +23.6 MeV, accounting for the Coulomb energy difference between the daughter nuclei 
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FIG. 3: Charge exchange SD strength dB<y ^ ^ (upper panel) and dB ^g + ^ (lower panel) of 90 Zr. 
The circles and squares are the experimental data taken from ref. 16]. The spectra dB ^^ + ^ are 
shifted by the Coulomb energy difference between the two daughter nuclei 90 Nb and 90 Y (+23.6 
MeV) to adjust the isospin difference between the two nuclei. The calculated results are plotted 
with the quenching factor quf = 0.68. The SD strength is averaged by the weighting function ([5]) 
with the width A = 2 MeV. 

90 Nb and 90 Yb. We introduce the quenching factor quf = 0.68 for both the t- and t + 
channels. For the t- channel, the experimental strength distribution peaked at E x ~ 26 
MeV is well described by the SLy4 interaction. The results of SGII and SIII also give 
reasonable agreement with the experimental peak energy. None of the calculated results 
show any substantial strength above E x ~ 36 MeV, while a significant portion of the sum 
rule value is found above E x ~ 36 MeV in the experimental data. This difference may be due 
to the lack of coupling to many-particle many-hole states in the present RPA calculations. 
In ref. 28|], the t_ SD strengths in 90 Zr have been studied using the RPA model including 
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the couplings to 2p-2h states. It was found that the mixing between lp-lh and 2p-2h states 
gives a large asymmetric spread in the strength of the SD resonances, and about 30% of the 
total strength is shifted to excitation energies above 35 MeV, referred to the parent nucleus 
90 Zr. This result is consistent with the quenching factor adopted in Fig. [3J It should be 
mentioned that the peak energy of the t_ SD strength is not changed appreciably by the 
coupling to the 2p-2h states, while the peak height is decreased substantially. 

For the t + channel, the two peak structures can be seen in both the calculated and 
experimental results. SkI3 and SLy4 describe the SD strength well at the low energy spectra. 
The calculated strength up to E x =40 MeV exhausts 100% of the sum rule value, while the 
experimental data show appreciable strength above E x =40 MeV. This difference may be 
due to the couplings to many-particle many-hole states similar to the t_ channel. 

Let us now discuss the integrated SD strength. The integrated SD strength 

ME,)- £ r d -§T dE ' ( 6 ) 



o 



At=0-,1-,2- 

is plotted as a function of the excitation energy E x in Fig. 0] for the operators = ^ t l _ 
7\\a x Y" 1 (f i )] A and Sw. = r A' y x ^i(j?)] A - The experimental data are taken from ref. 



161 ] . The value S_ is obtained by integrating up to E x = 50 MeV from the ground state 
of the daughter nucleus 90 Nb (E x = 57 MeV from the ground state of the parent nucleus 
90 Zr), while the corresponding value S + is evaluated up to E x = 26 MeV from the ground 
state of 90 Y (E x = 27.5 MeV from the ground state of the 90 Zr). This difference between 
the two maximum energies of the integrals stems from the isospin difference between the 
ground states of the daughter nuclei, i.e., T=4 in 90 Nb and T=6 in 90 Y. That is, the 23.6 MeV 
difference originates from the difference in excitation energy between the T=6 Gamow- Teller 
states in the (p,n) and (n,p) channels fl^. For both the S- and S + strength, the calculated 
results overshoot the experimental data in the energy range E x = 20-40 MeV. These results 
suggest the quenching of 30-40% of the calculated strength around the peak region, as was 
already mentioned. However, the integrated cross-sections up to E x = 56 MeV in Fig. H] 
approach the calculated values for both the £_ and t + channels. 

The calculated SD sum rule values in A=90 nuclei obtained by using the HF+RPA results 
are tabulated in Table HVl for the transitions with A 7r = 0~, 1~ and 2~. Clearly, the AS values 
show signs of multipole proportionality (2A+1), even though S"_ and S + themselves do not 
show any clear multipole dependence. The present RPA results for 90 Zr listed in Table IIVI 
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FIG. 4: Integrated charge exchange SD strength ([6]) excited by the operators S- = 
X)im pt l - a m r iYi and S + = m u t\v\n r iYi on 90 Zr. The calculated results are obtained 
by the HF+RPA model using the Skyrme interactions SIII, SGII, SLy4 and SkI3. The upper panel 
shows the S- and S+ strength, while the lower panel shows the S- — S+ strength. All strengths 
for the three multipoles A 7r =0 _ , and IT are summed up in the results. The experimental data 
are taken from ref. [lg]. No quenching factor is introduced in the calculation of the integrated 
strength. 

satisfy the sum rule value (2) in Table HT1 with high accuracy, to an error of only (0.1~0.2)%. 
This agreement guarantees the numerical accuracy of the present RPA calculations. This is 
also the case in 208 Pb, as will be shown in Section IIB. The AS = S- — S+ value is shown as 
a function of E x in the lower panel of Fig. HI We note that the AS* value saturates both in 
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TABLE IV: Sum rule values of charge exchange SD excitations in A=90 nuclei obtained by the 
HF+RPA calculations; S_ for 90 Nb and S+ for 90 Y. The SD strength is integrated up to E x = 50 
MeV for S_ and E x = 26 MeV for S+, respectively. The experimental data are taken from ref. 
lf|. The SD sum rules are given in units of fm 2 . See the text for details. 





SIII 


SGII 


SkI3 


SLy4 




S- 


S+ 


AS 


S- 


S+ 


AS 


S- 


s+ 


AS 


S- 


s+ 


AS 


0" 


34.8 


18.5 


16.4 


33.2 


17.4 


15.8 


36.6 


19.1 


17.5 


37.8 


21.4 


16.4 


1" 


120.8 


71.7 


49.1 


122.0 


74.3 


47.7 


120.8 


68.2 


52.7 


115.8 


66.4 


49.4 


2" 


130.1 


48.5 


81.6 


125.5 


45.9 


79.5 


139.0 


51.1 


87.9 


138.7 


56.4 


82.3 


sum 


285.7 


138.6 


147.1 


280.7 


137.6 


143.1 


296.3 


138.3 


158.1 


292.3 


144.2 


148.2 


exp 


S- = 271 ± 14 


S+ = 124 ± 11 


AS = 147 ± 13 





the calculated and the experimental values above E x = 40 MeV, while the empirical values 
S- and S+ themselves increase gradually above E x = 40 MeV. This is the crucial feature for 
extracting the model-independent sum rule AS* = S- — S + from the experimental data. The 
empirical values S + and AS" obtained from these analyses are shown in Table HVl The 
indicated uncertainties of SL, S + and AS contain not only the statistical error of the data, 
but also errors due to the various input of the DWIA calculations used in the MD analysis, 

n 

such as, the optical model parameters and the single-particle potentials [15|]. There is an 
additional uncertainty in the estimation of the SD unit cross-section, namely, the overall 
normalization factor 16|, which should be studied further experimentally. 

From AS, the neutron radius of 90 Zr is extracted to be \/< r 2 > n = (4.26±0.04)fm from 
the model- independent SD sum rule (j3J), where the empirical proton radius ^J<r 2 > p = 4.19 
fm is used. The proton radius is obtained from the charge radius in Table [Til by subtracting 
the proton finite size correction. The experimental uncertainty in the neutron skin thickness 



obtained by proton scattering is rather large: 5. 



up 



(0.09 ± 0.07)fm. This is 



because of the model-dependent analysis of the proton scattering, with effective nucleon- 
nucleon interactions in the nuclei [z]. On the other hand, the sum rule analysis of the 
SD strength determines the neutron radius with 1% accuracy, which is almost the same as 
that expected for the parity violation electron scattering experiment. The obtained value 
r n — f p = (0.07 ± 0.04) fm can be used to disentangle the neutron matter EOS by using the 
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aaa, 



strong linear correlation between the two quantities [l|, |2j, l3j , as will be discussed in Section 
3. 



B. Charge exchange SD excitations of Pb 



TABLE V: Proton, neutron and charge radii of 208 Pb. The charge radius is obtained by folding 
the proton finite size. The sum rule values AS = S- — S+ of the SD excitations are calculated 
by Eq. ([3]) with the HF neutron and proton mean square radii. Experimental data on the charge 



radius are taken from ref. 



scattering 



. Experimental data on 5, 



up 



and the giant dipole excitations of ^ u °Pb 



- r v are obtained by the proton 
The radii are given in units of 



fm, while the SD sum rules are given in units of fm 2 





SIII 


SGII 


SkI3 


SLy4 


exp 


r p 


5.521 


5.454 


5.421 


5.457 




r c 


5.578 


5.512 


5.479 


5.515 


5.503 ± 0.002 


r n 


5.646 


5.589 


5.649 


5.617 






0.125 


0.135 


0.228 


0.160 


0.083 < 5 np < 0.111[8], 0.19 ± 0.09[11] 


AS 


1086. 


1072. 


1154. 


1098. 





The HF results of 208 Pb are summarized in Table |V] The RPA results of SD excitations 
of 208 Pb are given in Figs. [5] and [6] for the four different Skyrme interactions, namely, SIII, 
SGII, SkI3, and SLy4. For the t_ channel, the strength distributions are spread out in a 
broad energy region (15 MeV< E x <35 MeV) except a tiny peak at E x ~ 5 MeV. On the 
other hand, the strength for the t + channel is concentrated in a single narrow peak. The 
highest peak of the t_ channel occurs at E x ~ 27-28 MeV in the cases of the SIII, SGII, and 
SLy4 interactions, while it is shifted to higher energies (E x ~ 33 MeV) in the case of SkI3. 
The 0~ and 1~ excitations have merged into one peak, having more than 40% of the total 
strength at the high energy side, while the 2~ states split into a broad energy region. The 
low-energy 2 _ state at around E x = 4 MeV is mainly due to the (vrl/ig^^l^^ 1 ^) excitation. 
The 0~ peak is predicted to occur at a slightly higher energy than the 1~ peak. However, 
it might be difficult to observe this peak experimentally because of its rather low strength. 
There are appreciable differences in the peak energies between the Skyrme interactions for 
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FIG. 5: (Color online) Charge exchange SD strengths for the operators = r ii a x 

calculated by the HF+RPA model with the Skyrme interactions (a) SIII, (b) SGII, (c) SkI3, and 
(d) SLy4. The excitation energy is referred to the ground state of the parent nucleus 208 Pb. The 
SD strength is averaged by the weighting function in Eq. ([5]) with the width A = 1 MeV. 



the t+ channel: E x ~ 3 MeV for SGII, E x ~ 5 MeV for SIII and E x ~ 6 MeV for SLy4 and 
SkI3, as listed in Table IVO The sum rule values S_ and S + are listed in Table [VTT1 Because 
of the strong Pauli blocking of neutron excess in 208 Pb, the S + value is much smaller than 
the S_ value, at most, 20% of the corresponding 5_ value for each multipole. The S+ value 
is substantial in the case of A=90 as shown in Table I IV} more than 55% of S_ in some cases. 
However, AS = S_ — S + obeys the (2A+1) proportionality, as expected from Eq. (j2J). The 
charge exchange 208 Pb( 3 He,t) 208 Bi reaction was performed to study the SD strength in 208 Bi. 
The data were analyzed by a least-squares fitting method and the peak of the SD strength 
was found to be at E x = 24.8±0.8 MeV, as measured from the ground state of 208 Pb |29| . 
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FIG. 6: (Color online) Charge exchange SD strengths for the operators = r i[ a x 

calculated by the HF+RPA model with the Skyrme interactions (a) SIII, (b) SGII, (c) SkI3 and 
(d) SLy4. The excitation energy is referred to the ground state of the parent nucleus 208 Pb. The 
SD strength is averaged by the weighting function in Eq. ([5]) with the width A = 1 MeV. 

This empirical peak energy is close to the average energy E of SD strength obtained by 
SIII and SGII in Table IVII Further experimental effort is urgently needed to obtain more 
quantitative strength distributions, for example, for the multipole decomposition analysis of 
charge exchange reactions on a 208 Pb target. 

One can see only one sharp peak in the t + channel in Fig. [61 There are only two allowed 
lp-lh configurations (^gg^nlh^^) and (ulin/27Tlh^y 2 ) for both 1~ and 2~ excitations 
because of the strong Pauli blocking effect of excess neutrons. Moreover, the v2g^/2 and 
z/li n /2 states are almost degenerate in energy in the HF potential. They are the reasons why 
there is only one sharp peak in the t + channel of 208 Pb. It might be interesting to perform 
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FIG. 7: Charge exchange SD strength dB<y ^~^ (upper panel) and dB ^^ + ^ (lower panel) of 208 Pb. 
The spectra dB<y ^ are shifted by +37.2 MeV due to the Coulomb energy difference between the 
two daughter nuclei 208 Bi and 208 T1. The arrow in the upper panel shows a peak energy at E x = 
24.8 MeV observed by the charge exchange reaction 208 Pb( 3 He,i) 208 Bi [29]. 



208 Pb(n,p) 208 Tl or 208 Pb(t, 3 He) 208 Tl reactions in order to observe this peak experimentally. 
The 208 Pb(n,p) 208 T1 reaction has been reported for the t + channel, and a broad peak found 



at E x ~ 8MeV, as measured from the ground state of 208 Pb with rather poor statistics 30]. 

The integrated SD strengths for both the £_ and t + channels are shown in Fig. [HJ The 
calculated NEWSR shows a saturation at around E x ~ 30 MeV as can be seen in Fig. [HJ 
As noted previously, the t + channel has only a small contribution to the model-independent 
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FIG. 8: Integrated charge exchange SD strength (J6j) of 208 Pb for the operators S- = 
^i,m,ut- a m r iYi(ri) and S + = Yli,m,u^+ a m r i^'l (^i) calculated by the HF+RPA model with the 
Skyrme interactions SIII, SGII, SkI3 and SLy4. The upper panel shows the S- and S+ strength, 
while the lower panel shows the AS = S- — 5+ strength. All strengths for the three multipoles 
A 7r =0 _ , 1~ and 2~ are summed up in the results. 

sum rule AS. 

The couplings to the 2p-2h states may increase the spread in the SD strength in A=208 
nuclei as well as A=90 nuclei. So far, the charge exchange Gamow-Teller(GT) states in 208 Bi 
were studied by taking into account the couplings to 2p-2h states in the particle-vibration 



model 



24| . While a large spread was found in the GT states in the particle- vibration model 



18 



TABLE VI: Peak energies and the average energies of charge exchange SD excitations in A=208 
nuclei calculated by the HF+RPA model; S_ for 208 Bi and S+ for 208 T1. The average energy is 
calculated by the ratio of EWSR to NEWSR: E(MeV)=mi/m . See the text for details 





t- 


t+ 




E peak {MeV) 


E(MeV) 


E peak (MeV) 


E(MeV) 


SIII 


26.7 


24.2 


5.0 


7.3 


SGII 


28.1 


24.6 


2.5 


6.0 


SkI3 


32.7 


27.9 


5.6 


7.3 


SLy4 


27.4 


23.6 


6.3 


8.0 



TABLE VII: Sum rule values of charge exchange SD excitations in A=208 nuclei calculated by the 
HF+RPA model; S_ for 208 Bi and S+ for 208 T1. The SD strength is integrated up to E x = 57 MeV 
for S_ and E x = 20 MeV for S+; the excitation energy is referred to the ground state of 208 Pb. 
The SD sum rules are given in units of fm 2 . See the text for details. 





SIII 


SGII 
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S- 


s+ 


AS 


S- 


s+ 


AS 


S- 


s+ 


AS 


S- 


s+ 


AS* 


0" 


148.6 


27.0 


121.6 


114.1 


24.3 


119.8 


158.0 


29.7 


128.3 


158.5 


36.0 


122.5 


1- 


442.7 


78.8 


363.9 


440.4 


82.3 


358.1 


454.5 


69.2 


385.3 


430.8 


63.6 


367.2 


2" 


632.2 


28.3 


603.9 


620.7 


26.4 


595.3 


669.8 


28.2 


641.6 


644.5 


34.1 


610.5 


sum 


1224. 


134.1 


1089. 


1205. 


132.0 


1073. 


1282. 


127.1 


1155. 


1234. 


133.7 


1100. 



calculations, the peak energy did not change appreciably due to the couplings to 2p-2h 
states. There have been no microscopic studies of SD states that take into account the 
couplings to 2p-2h states in A=208 nuclei. 



III. SD SUM RULES AND NEUTRON MATTER EOS 

Sum rules are useful tools to study the collective nature of excitation modes in many- 
body systems. In particular, for charge exchange excitations, model-independent sum rules 
are derived and used to analyze experimental data on Gamow- Teller resonances and SD res- 
onances [yl . For SD states, the sum rules can be used to extract the neutron skin thickness, 
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as was discussed in Section 2. References 

m 

have reported a strong correlation between 
the neutron skin thickness and the neutron matter EOS, as obtained by using Skyrme and 
relativistic mean field theories. In this section, we will study the relation between the SD 
sum rules and the neutron matter EOS. The strong linear correlation between the neutron 
skin thickness 

S nP = VWn - (7) 

and the pressure of neutron matter 

d(E(p n )/p n ) 

dpn 

is essential for this study. Other linear correlations between the neutron skin thickness 
and various isovector nuclear matter properties have also been pointed out recently {(J. 
Given these correlations, accurate information on the neutron skin thickness will be quite 
useful in determining empirically the pressure of neutron matter EOS and isovector nuclear 
properties, such as the volume and surface symmetry energies. 

The correlations between the pressure of neutron matter at the neutron density p n = 0.1 
fm -3 and the charge exchange SD sum rules of 90 Zr and 208 Pb are shown in Figs. [9] and 
[10] with 12 different Skyrme interactions. The numbers denote different Skyrme parameter 
sets: 1 for SI, 2 for SIII, 3 for SIV, 4 for SVI, 5 for Skya, 6 for SkM, 7 for SkM*, 8 for 
SLy4, 9 for MSkA, 10 for SkI3, 11 for SkX and 12 for SGII. The correlation coefficients 
from the extrapolated lines are r = 0.888 and 0.811 for 208 Pb and 90 Zr, respectively. The 
correlation coefficients are somewhat smaller than those of the calculated correlation between 
the neutron skin thickness 5 np and the pressure P in ref. but still, we can see fairly good 
correlations in Figs. l9land[T0l 

The rms proton, charge, and neutron radii in 90 Zr calculated by the HF model with the 
four interactions SIII, SGII, SkI3 and SLy4 are shown in Table HJJ The calculated charge 
radii of the SGII and SkI3 interactions show reasonable agreement with the experimental 
values. However, there is a factor 2 difference in the neutron skin thickness 5 np between the 
two interactions. As seen in Table [TTJ the neutron skin thickness 5 np obtained by the SD 
sum rules is consistent with the value previously obtained from the proton scattering data. 
However, the experimental uncertainty in the value S np = (0.07 ± 0.04) fm obtained by the 
SD sum rules is half that obtained through the proton data. This small uncertainty will help 
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FIG. 9: Correlations between the pressure of neutron matter and the SD sum rule values of 90 Zr 
with 12 different Skyrme interactions. The numbers denote different Skyrme parameter sets: 1 for 
SI, 2 for SIII, 3 for SIV, 4 for SVI, 5 for Skya, 6 for SkM, 7 for SkM*, 8 for SLy4, 9 for MSkA, 10 
for SkI3, 11 for SkX and 12 for SGII. The correlation coefficient is found to be r = 0.811. 



to disentangle the neutron matter EOS using the strong correlation with the neutron skin 
thickness. The experimental skin thickness 5 np = 0.07 ± 0.04fm is close to the HF results of 
SLy4, as well as SGII and SIII. The SkI3 result is not favored over the empirical result, even 
taking the experimental uncertainties into consideration. We should also note that the ex- 
perimental peak energy of t_ SD strength in 90 Nb coincides with the calculated peak energy 
of the SLy4 interaction, while that of SkI3 is 4 MeV above the experimental value, as seen 
in Fig. [3] While all interactions lie within the experimental value AS = (147 ± 13)fm 2 in 
Fig. [9l the empirical data favor the interactions indicated by the numbers 2(SIII), 11 (SkX), 
8(SLy4), 7(SkM*) and 6(SkM). These interactions suggest a soft neutron matter EOS with 
the pressure P(p„=0.1fm~ 3 ) = (0.65±0.2) MeV. Thus, the preferred nuclear matter symme- 
try energy extracted from the SD experiment is found to be J = (30±2) MeV as a result of 
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FIG. 10: Correlations between the pressure of neutron matter and the SD sum rule values of 208 Pb 
with 12 different Skyrme interactions. The numbers denote different Skyrme parameter sets: 1 for 
SI, 2 for SHI, 3 for SIV, 4 for SVI, 5 for Skya, 6 for SkM, 7 for SkM*, 8 for SLy4, 9 for MSkA, 
10 for SkI3, 11 for SkX and 12 for SGII. The dashed line represents the result obtained by the 
least-squares method. The correlation coefficient is found to be r = 0.888. 

the strong correlation between the neutron skin thickness and the symmetry energy 

Table |V] tabulates the rms proton, charge and neutron radii in 208 Pb calculated by the 
HF model, along with the experimental charge radius. The HF results of SGII and SLy4 
account for the experimental charge radius, while there is a large variation in the predictions 
for the neutron skin thickness 5 np . The empirical value of the neutron skin thickness 5 np 
in 208 Pb was obtained by proton scattering experiments. However, the values obtained 
depend very much on the experiments and analyses. That is, the experimental errors are 
still large and some of the values obtained have no overlap, even when the uncertainty in 
the analyses is taken into account; 5 np = (0.14 ± 0.02) fm in ref. 5 np = (0.20 ± 0.04) 
fm in ref. Q and (0.083 < S np < 0.111) fm in ref. [8|. We quote in Table El the value 
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in ref. Q where the ana>yses were performed eo m prehe„ 8 ive. y with many different 8 ets of 
data including those adopted in refs. |9J, llQj. Although these results depend on the effective 
nucleon-nucleon effective interactions in nuclei used in the analysis, the comprehensive study 
of proton scattering in ref. [§] reports rather small neutron skin thicknesses, even smaller 
than the smallest value in Table IVl obtained using the SIII interaction. Again, this small S np 
suggests a soft neutron matter EOS similar to the conclusion reached by the SD sum rules 
of 90 Zr. The charge exchange 208 Pb( 3 He,t) 208 Bi reaction data [29j show an SD peak in 208 Bi 
at E x = 24.8±0.8 MeV measured from the ground state of 208 Pb, as marked by an arrow 
in Fig. [7J This peak position is close to the calculated value of the SGII interaction, while 
the SkI3 peak is a few MeV higher than the empirical value. This comparison may exclude 
the prediction by SkI3, which gives a hard neutron matter EOS in Fig. [TU] marked by the 
number 10. 

The neutron skin thickness was determined by the giant dipole resonance experiment 
to be 5 np = (0.19 ± 0.09)fm [111] . This analysis depends on the adopted transition density 
and also the optical potentials so that the result is highly model-dependent. We definitely 
need more quantitative information, i.e., model-independent information on the neutron 
skin thickness in 208 Pb for precise determination of the neutron matter EOS as well as the 
isovector nuclear matter properties. To this end, the charge exchange SD experiments of 
208 Pb will provide useful model-independent information with the same accuracy as the 
parity violation electron scattering experiment. 



IV. SUMMARY 

We have investigated the SD excitations in 90 Zr and 208 Pb using the HF + RPA model 
with four Skyrme interactions, viz., SIII, SGII, SkI3 and SLy4. It is shown that the Lan- 
dau damping effect plays an important role in explaining the large observed width of SD 
resonance, while the coupling to the continuum is rather weak. Among the four interac- 
tions, the peak position of the experimental £_ SD strength in 90 Nb is well described by 
the SLy4 interaction, while the results of SIII and SGII are also acceptable. For the t + 
excitation of 90 Zr, a two-peak structure was found in both the experimental and calculated 
results. The SLy4 and SkI3 results showed good agreement with the observed low energy 
peak. We pointed out that the calculated results need a quenching factor quf~0.68 to al- 
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low a quantitative comparison with the experimental data up to E x = 36(40) MeV for the 
£-(£+) channel in Fig. [31 About 30% of the NEWSR value is found in the excitation energy 
above E x =36(40) MeV for the 90 Zr(p,n) 90 Nb ( 90 Zr(n,p) 90 Y) experiments. The calcu- 
lated SD sum rule AS = 5_ — S + shows good saturation properties above E x =40 MeV 
without any quenching factor relative to the observed data despite the fact that sum rules 
S- and S+ themselves increase gradually above E x >40 MeV. The neutron skin thickness 
&np — 0.07 ± 0.04 fm extracted from the SD sum rules is close to the calculated values ob- 
tained using SLy4 as well as SIII and SGII. However, the extracted value does not favor the 
SkI3 interaction which gives almost twice as large a neutron skin thickness as SIII and SGII. 
This is indicative of the soft neutron matter EOS induced by the strong linear correlation 
between the neutron matter EOS and the neutron skin thickness. We showed that the SD 
strength of the t_ excitation of 208 Pb has a large width due to the Landau damping effect. 
In contrast, the t + excitation of 208 Pb turns out to be a single peak in a rather low energy 
region because of the strong Pauli blocking effect of the excess neutrons. The peak of the t- 
SD strength was observed by 208 Pb( 3 He,t) 208 Bi at E x ~25 MeV. This peak energy coincides 
with the peak calculated using the SGII interaction, while the SkI3 interaction yields a peak 
that is a few MeV higher than the empirical peak. Thus, the empirical SD sum rule values 
of 90 Zr and the observed peak energies of the i_ SD strength distributions in 90 Nb and 208 Bi 
indicate a soft neutron matter EOS with a pressure of P(p n =0.1fm -3 ) = (0.65±0.2) MeV. 
The nuclear matter symmetry energy is also determined to be J = (30±2) MeV from the 
strong correlation between the neutron skin thickness and the symmetry energy. In order to 
draw a more definite conclusion on the SD sum rules, as well as the neutron skin thickness 
and the neutron matter EOS, we need quantitative experimental work to obtain the SD sum 
rules in heavy nuclei like 208 Pb, both in the t- and t + channels. 
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